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^'- ABSTRACT 

We present a comprehensive survey of C 11* absorption detections toward stars 
within 100 pc in order to measure the distribution of electron densities present in 
f^ I the local interstellar medium (LISM). Using high spectral resolution observations 

of nearby stars obtained by the Goddard High- Resolution Spectrograph (GHRS) 
and the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space 



O 



c^ I Telescope (HST), we searched for all detections of LISM C 11* absorption. We 

identify 13 sight lines with 23 individual C 11* absorption components, which 
_ provide electron density measurements, the vast majority of which are new. We 

CNj I employ several strategies to determine more accurate C II column densities from 

^~ I the saturated C II resonance line, including, constraints of the line width from the 

optically thin C II* line, constraints from independent temperature measurements 
"^ I of the LISM gas based on line widths of other ions, and third, using measured S II 

lyQ I column densities as a proxy for C II column densities. The distribution of electron 

^ I densities based on using S II as a proxy for C II is similar to the distribution 

based on carbon alone, while significantly tighter, and proves to be a promising 
/\ ■ technique to avoid grossly overestimating the C II column density based on the 

c3 ] saturated line profile. The sample of electron densities appears consistent with a 

log- normal distribution and an unweighted mean value of ne{C IIsii) = 0.11!^o;q5 
cm~^. Seven individual sight lines probe the Local Interstellar Cloud (LIC), 
and all present a similar value for the electron density, with a weighted mean of 
ne(LIC) = 0.12±0.04 cm~^. Two clouds, the NGP and Gem clouds, show similar 
electron density properties as the LIC. The Hyades Cloud, a decelerated cloud at 
the leading edge of the platoon of LISM clouds, has a significantly higher electron 
density than the LIC. Observed toward G191-B2B, the high electron density 
may be caused by the lack of shielding from such a strong radiation source. 
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No evidence of a correlation between electron density and angular separation of 
the sight line from the strongest extreme ultraviolet radiation source, e CMa, is 
found. Given some simple assumptions, the range of observed electron densities 
translates into a range of thermal pressures, P/k = SSOOI^qqq K cm~^. This 
work greatly expands the number of electron density measurements and provides 
important constraints on the ionization, abundance, and evolutionary models of 
the local interstellar medium. 



Subject headings: atomic processes — ISM: abundances — line: profiles — solar 
neighborhood — techniques: spectroscopic — ultraviolet: ISM 



Introduction 



More than a decade of use of high spectral resolution {R = X/AX > 50,000) ultra- 
violet (UV) spectrographs with wide spectral coverage, such as the Goddard High Reso- 
lution Spectrograph (GHRS) and particularly the Space Telescope Imaging Spectrograph 
(STIS) onboard the Hubble Space Telescope (HST), has been a boon for observations of the 
local interstellar medium (LISM). Since the LISM, the collection of warm gas in the im- 
mediate (<100 pc) vicinity of the Sun, is traversed by all sight lines that extend beyond 
our solar system, the number of observations with LISM absorption has grown tremen- 
dously. The use of this expanded database has been the impetus for many investigations 
into the nature of the LISM, including the other papers in this series, which present fits 
to LISM absorption in many i ons (e.g., D I, C II, N I, O I, Mg II, Al II, Si II, and Fe II, 
Redfield fc Linskvll2002l. l2004al). t emperature and turbulent velocity mea surements in LISM 
cloud s (JRedfield fc Linskyll2004bl ). and a dynamical model of the LISM (JRedfield fc Linsky 
2008h . 



Due to the low column densities (logA^(H I) ~ 16.8 — 18.3) typical of local clouds, 
the ionization structure of hydrogen is vital to understanding the physical structure and 
origins of the LISM. An accurate accounting of the ionized fraction of hydrogen (and the 
ionization levels of all ions) is critical for measuring abundances and the depletion of gas 
phase ions onto dust grains. Low column density clouds, (i.e., logA^fH I) < 19.5), like the 



LISM clouds, are not signific antly shielded to ionizing photons (j Jenkins! |2004J : iJenkins et al. 



^1 -n^-\ ^ - ^ 

2000bl : iSofia fc Jenkins! 119981 ). The origin and evoluti on of the Local Bubbl e, the ~100 pc 

radius cavity in which the warm LISM clouds reside (JLallement et al.ll2003l ). is encoded in 

the ionization structure. Recent nondetections of high temperature lines in the extreme 



UV flHurwitz et a. 



sphere (JLallement 



20051) and the realization that soft X-ray emission caused by the helio- 



2004! ) may contribute to the emission formerly assigned to nearby hot gas 
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( jSnowden et al.lll990l ). highlight the current challenges in understanding the thermal struc- 
ture and ionization level of the Local Bubble. A realistic Local Bubble can be modeled far 



from ionization equilibrium (JBreitschwerdt fc de Avillezll2006l ) and the physical structure of 
the Local Bubble can have a direct influence on the ionization structure of the warm LISM 



clouds (ISlavin fc FrischI 120021 ) . Additionally, the structure of the heliosphere, the interface 
between the LISM and the solar win d, depends signific antly on the ionization structure of 
the surrounding interstellar medium (JMiiller et al.ll2006l ). 



The ionization fraction of hydrogen, X{B.) = n(H II)/(r2(H I) + n(H II)), is computed 
using the ?t,(H I) = n(He I) x A^(H I)/A^(He I) and n(H II) ^ Ug. Typically, ?7,(He I) is derived 
from in situ measurements o f interstellar helium st reaming through the heliospheric interface 
into the inner solar system (JGloeckler et al.ll2004l ). which of course, is only a measurement 
of the helium density at one point in the LISM. A^(H I)/A^(He I) is derived from extreme 
UV observations of nearby white dwarfs (WDs), which contain the ionization edge of He I, 
and the continuum provides an estimate of H I (JDupuis et al.l Il995l : iBarstow et al.l 119971 ). 
The weighted mean from this work, based on nine white dwarf sight lines within 100 pc, is 
A^(H I)/iV(He I) = 14.1 ± 1.7, although this is a full sight line average, since this technique 
cannot separate individual cloud components. 

The remaining measurement in the calculation of the ionization structure is the elec- 
tron density, rig. The methods employed to make this measurement utilize atomic transitions 
along interstellar sight lines, which have the benefit of being able to resolve individual ab- 
sorbers if taken at high spectral resolution and can provide a large number of measurements 
through various LISM environments. The ratio of magnesium ionization stage column den 



sities , A^fMg ll)/N(M^ I), has provided a number of rip, rneasur ements (e.g., iFrisch et al. 



1990l : iLallement et al.lll994j : lFrischlll994j : iLallement fc Ferletlll997l ). However, this technique 
suffers from the requirement of ionization equilibrium and a strong temperature dependence. 
Alternatively, the ratio of the collisionally excited carbon line column density to the reso- 
nance line column density, A^(C II*)/A^(C II), can provide rig estimates without the need for 
ionization equilibrium and has a very weak temperature dependence. However, due to the 
weakness of the excited C II* absorption and saturation of the only availab le UV C II reso- 
nance line, few LISM sight lines have been analyzed using this technique. IWood &: Linsky 
(119971 ) used high spectral resolution {R ~ 100,000) spectra of a Aur to calculate He along 
the line of sight and demonstrated that the relatively simple absorption profiles through the 
LIS M provide an excellen t opportunity to make precise measurements of the electron den- 
sity. iHolberg et al.l (119991 ) used the same technique to measure the electron density in ISM 
absorbers along the line of sight toward WD 1029+537, a white dwarf ~132pc away. The 
excited absorption along other LISM sight lines has been ra easured at moderate resolution 



{R ~ 20,000) with Copernicus (e.g., I York fc Kinahanlll979l ) and the Far Ultraviolet Spec- 
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troscopic Explorer {FUSE) (JLehner et al.ll2003l ). although typically individual absorbers are 
not resolved, so the estimates are full sight line averages. A few rig measurements in slightly 
warmer, and possi bly local, O Vl-b earing gas, derived by comparisons of O VI emission 
(JDixon et al.l l2006l ) and absorption (JLehner et al.ll2003l ). have led to rif. estimates of ~0.2 



cm 



We present an inventory of high spectral resolution C II* detections along LISM sight 
lines from the complete HST spectroscopic database. We employ various strategies to cir- 
cumvent one of the challenges of using the C ir/C II ratio technique, namely obtaining a 
reliable C II column density from a saturated resonance line: (a) Simultaneously fitting the 
optically thin C 11* and saturated C II profiles puts realistic constraints on the Doppler width 
and column densities of the resonance line; (b) Use of previously determi ned LISM temper- 
atures and turbulent velocities, typically derived from the same dataset (JRedfield &: Linsky 
2004bl ). constrains the Doppler width of the carbon lines; and (c) Use of optically thin pro- 
files of S II, which has a sirnilar i onization potential as C II, as a proxy of the C II column 
density (e.g., lOliveira et al.l 120031 ). This collection of new, high spectral resolution UV ob- 
servations of nearby stars presents an opportunity to greatly expand the number of LISM 
electron density measurements and probe the ionization structure of our most immediate 
interstellar environment. 



2. Observations 



For our purpose, we are interested in sight lines toward stars within 100 pc that show 
interstellar absorption in both C II (1334.5323 A) and the C 11* doublet (1335.6627 A and 
1335.7077 A). We compiled all moderate to high resolution observations of nearby stars with 
the i/ST spectrographs: GHRS and STIS. The complete sample includes 417 unique targets 
within 100 pc, almost half of which have spectra that cover the wavelength region of C II 
and C 11*. All these relevant spectra were scrutinized for signs of C 11* absorption. We found 
only 13 sight lines that show LISM absorption in both transitions and list them in Table [H 

The 100 pc distance limit is chosen to coincide with the approximate extent of the Local 
Bubble ( jLallement et al.ll2003l ). Observations of more distant stars may be more difficult to 
analyze, as they are more likely to traverse many absorbing clouds creating a blended line 
proffie. 

Table [2] lists the observational parameters of the datasets extracted from the HST'DaXa, 
Archive. All observations are necessarily moderate {R > 20,000) to high [R > 100,000) 
spectral resolution in order to resolve narrow, closely spaced interstellar absorption features 
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and to increase the likelihood of the detection of the weak C 11* line. 

We analyzed LISM absorption in the S II multiplet (1250.578 A, 1253.805 A, and 
1259.518 A) in order to better constrain the C II column density. Ten of the 13 targets 
have spectra that include the S II lines. In addition, we fit the LISM absorption in Mg I 
(2852.9631 A) and the Mg II doublet (2796.3543 A and 2803.5315 A) in order to further 
constrain the electron density. Only three of the 13 targets have spectra that include both 
Mg I and Mg II. Observations utilized for these analyses are also included in Table [21 

We reduced the GHRS data acquired from the HST Data Archive wi th the CAL HRS 



software package using the Image Reduction and Analysis Facility (IRAF; lTodylll993l ) and 
the Space Telescope Science Data Analysis System (STSDAS). We used the most recent 
reference calibration files. Many of the echelle observations were obtained in the FP-SPLIT 
mode to reduce fixed-pattern noise. The individual readouts of the FP-SPLIT spectra ar e 



combined using a cross-correlation procedure called HRS_MERGE (JRobinson et al.lll992l ). 
The reduction included assignment of wavelengths using calibration spectra obtained during 
the course of the observations. The calibration spectra include either a WAVECAL, a direct 
Pt-Ne lamp spectrum used to derive the dispersion relation, or a SPYBAL (Spectrum Y- 
Balance), which only provides a zero-point wavelength offset. Any significant errors involved 
in the wavelength calibration are included in our central velocity determinations. The wave- 
length calibration of the Fe II spectrum of a Gru is used to calibrate the wavelength solution 
of C II, since no wavelength calibration of this segment of the spectrum was taken at the 
time. 

We reduced the STIS data acquired fro m the HST Da ta Archive using the STIS team's 



CALSTIS software package written in IDL (JLindlerlll999l ). The reduction included assign 



ment of wavelengths using calibration spectra obtained during the course of the observations. 
We used the ECHELLE_SCAT routine in the CALSTIS software package to remove scat- 
tered light. However, the scattered light contribution is negligible in this spectral range, and 
does not influence the uncertainties in our spectral analysis. 



3. Data Analysis and Line Profile Fitting 

Figures [la] and [lb] show the C II resonance and the C 11* excited absorption lines 
observed toward the 13 targets listed in Table [U The spectra (histogram) are plotted in 
heliocentric velocity. Also plotted are the individual interstellar component fits (dashed 
lines) and the total interstellar absorption convolved with the instrumental profile (thick 
solid lines). As explained in Section [3^21 the resonance line plots come from the simultaneous 
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fits, and the excited line plots come from the individual excited line fits. 

The stellar continuum level (thin solid lines) illustrates our flux estimates of the spectra 
without any interstellar absorption. These continua were determined by fitting polynomials 
to the spectral regions adjacent to the observed interstellar features. If the intrinsic spectrum 
is a complex profile with intricate stellar features, placing the continuum can be a difficult 
task. All of our target spectra have relatively smooth continua, and we are able to easily 
reproduce the intrinsic continuum flux backgrounds with low order polynomials. Only two 
of the objects analyzed (a Aur and EX Hya) include late-type stars, which are characterized 
by stellar emission line spectra as opposed to smooth continuum spectra (e.g., as seen in 
the white dwarf spectra), although their broad emission lines allowed for simple continuum 
placement. 

We used a Gaussian ISM absorption profile fitting algorithm to de termine t he fits . The 



program makes use of rest wavelengths and oscillator strengths from iMortoru (120031). and 



the in stru mental line s p read functions for GHRS and STIS spectra are taken from iGilliland 



(119941 ) and lSahuet all (I19991 ). respectively. 



3.1. Fit Parameters 

Table[3]lists the fit parameters measured, including central velocity {v [km s^^]), Doppler 
width [h [km s~^]) and the logarithm of the column density (A^ [cm~^]). The central velocity 
is the mean radial velocity of the absorbing component. The Doppler width is a function of 
the temperature (T [K]) and turbulent velocity (^ [km s^l]) of the interstellar gas: 

52 = f^ + ^2^ 0.016629^ + e', (1) 

m A 

where k is Boltzmann's constant, m is the mass of the observed ion, and A is the atomic 
mass {Aq = 12.011). The column density is a measure of the amount of material along the 
line of sight to the target. 



3.2. Procedure 

For each target, we first fit the excited line alone. Since the excited line is actually a 
doublet, each paired absorption line possesses the same intrinsic central velocity, Doppler 
width, and column density; so these parameters are constrained accordingly to be identical 
for both lines of the doublet. We then fit the excited (C 11*) and resonance (C II) lines 
simultaneously, requiring the central velocities and Doppler widths to be identical. Since 
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the resonance line is typically saturated, A^(C II) is not well constrained, but fitting the 
lines simultaneously allows us to make use of information retained in the, albeit saturated, 
resonance line, while retrieving vital constraints from the optically thin excited line. Thus, 
the final fit parameters for the resonance line are derived from the simultaneous fits only. 
For the excited line fit parameters, we supplement the results from the simultaneous fits 
with those from the individual fits in order to limit the impact of systematic errors; these 
values are weighted means of the parameters from both the simultaneous fits and individual 
excited line fits. 

In saturated lines, Doppler width and column density are tightly coupled. Thus, any 
available constraints on the b value will yield a more certain measurement for N, which in 
the case of C II and C 11*, is critical in determining electron d ensities. Consequently, we use 
measured LISM temperatures and turbulent velocities from iRedfield &: LinskyI (l2004bl ) to 
constrain the Doppler widths for the a Aur, rj UMa, and G191-B2B lines of sight. Unfortu- 
nately, LISM temperatures and turbulent velocities have been determined for only a limited 
number of sight lines, so we are able to include this type of constraint for only these three 
targets. 

Figure [lb] shows four white dwarfs (WD 0050-332, WD 0232+035, WD 2309+105, and 
WD 1210+533) that each possess a resonance line component that does not have an excited 
line counterpart. In each case, the Doppler width and column density for the lone resonance 
component are largely determined by what is not accounted for by the other paired com- 
ponents; the C II measurements of the component without a C 11* counterpart are poorly 
determined as they lack the critical constraints provided by the optically thin excited com- 
ponent. Upper limits (3cr) are determined for the C 11* column density for these components 
that show absorption only in the resonance line. 

The dominant source of systematic error, as mentioned earlier, is the saturation of the 
C II reson ance line. It is made evident upon comparing o ur results for WD 232+035 with 
those from lVennes et al.l (120001 ). Using the same datasets, IVennes et al.l (l2000l ) observed two 
interstellar components toward WD 0232+035 with central velocities of 3.1 ± 0.2 km s~^ 
and 17.6 ± 0.9 km s"^, in good agreement with our measurements of 3.81 ± 0.47 km s~^ 
and 17.4 ± 4.6 km s~^. Their measured column density of logA^(C II) = 14.32 ± 0.57 
is a factor of 7 lower than our measurement of logA^(C II) = 15.18 ± 0.42, but both have 
large error bars, and the difference is just over Icr. In contrast, the resonance line comp onent 
that has no accompanying excited line component is measured by lVennes et al.l (120001 ) with 
a column density logA^(C II) = 14.16 ± 0.41, a factor of 60 lower than our measured 
value of logiV(C II) 



15.961q;55, ^^d a difference of almost 3a. The critical factor in the 



disagreement concerns the treatment of the Doppler width (6). IVennes et al.l (j2000l ) allow for 



a range of b values (3.5-10.0 km s~^ and 4.0-10.0 km s^^ for the ~3 km s^^ and ~17 km s"^ 
components, respectively), while we allow b to remain an independent variable. Such a wide 
range of b values appears to be unnecessary; for our C II fits to all sight lines, the weighted 
mean Doppler width is 6 = 4.30 ± 0.87 km s~^, and inconsistent with a large b value for 
WD0232+035. In addition, our fit to the S II lines, which also shows two components at 
similar velocities (4.11 ± 0.64 km s~^ and 16.8 ± 3.8 km s~^), has a weighted mean Doppler 
width of 6 = 3.95 ± 0.79 km s~^, also inconsistent with a large Doppler width. 

Simultaneously fitting the unsaturated excited fine with the accompanying saturated 
resonance line is the critical property of our method, which allows us to more tightly constrain 
the observed column density of the saturated line than if we were fitting it alone. Essentially, 
the optically thin excited line provides the constraint on the Doppler width, which thereby 
constrains the acceptable range of column densities. 



3.3. Component Determination 

In determining the number of absorption components, we use the fewest number of 
components that produce a satisfactory fit. In most cases, the resolution is high enough 
and the differences between component velocities are great enough such that discernment of 
components is almost trivial. Other times, however, it is not so easily determined. Com- 
ponents are added if a clear asymmetry is detected in either the C II and C 11* profiles 
(e.g., 7] UMa, WD 1210+533, IX Vel), or if other ions indicate that additional components 
are required (e.g., a Gru, p Lup). Our general attitude toward this part of the analysis is 
to approach the absorption feature with Ockham's razor in hand, so as to cut away any 
unjustified components. 



3.4. S II 

We are motivated to find an alternative means of estimating the C II column density, 
due to the difficulty of obtaining it directly from the strongly saturated resonance line. We 
use the optically thin S II triplet, located near the C II lines, as a proxy for C II due 
to their similarity in ionization potential, which is 24.4 eV for C II and 23.3 eV for S II, 
an d their similarity in th e ratio of their ionization and recombination rates (defined as P 



m 



Sofia fc JenkinsI (Il998l ). where S II and C II, among other ions, are compared in their 
Figure 3). Ten of the 13 sight lines with C 11* absorption also have spectra that cover the 
S II wavelength range. We fit all three S II lines simultaneously, and due to the range of 
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opacities of the three optically thin lines, we are able to measure an accurate S II column 
density. The spectra and fits are shown in Figures [2a] and [2bl and the fit parameters are 
listed in Table HI 

In order to convert from A^(S II) to A^(C II), we need to take into account the different 
abundances and depletion levels of sulfur and carbon in the LISM: 



A^(CIIsii) = N{Sll) X io[C0+^{c)]-[S0- 



-D(S)] 



(2) 



where A^(C IIsii) is the estimated C II column density based on S II as a proxy, A^(S II) is 
our measured column density of S II, C,^ = 8.39 ± 0. 05 and Sq = 7.14 ± 0.05 are the solar 



abund ances of carbo n and sulfur (JAsplund et al.ll2005l ). and D{C) and -D(S) are the depletion 



l evels (iJenkins 



( iLehner et al 



2004 ) . Here Z)(S)~ 0.0 , as has been detected in local and more distant ISM 
2OO3I : IWeltvet al.lll999h . and D{C) 



-0.17 ± 0.19 from iJenkinsI fl2004[ ). 
where the error inc orpora tes some of the na tural v ariation in the ISM (e.g., -D(C) ~ 0.0 



CardeUi et al.lll996l : — 0.2 lLehner et al.ll2003l : — 0.4 lWeltvet al.lll999f ). 



The C II column density estimated using this technique is listed in Table [3], and the 
comparison with the direct C II column density based on the C II and C 11* fits shown in 
Figures [Ta| and [lb] is exhibited in Figure [1] We present two errors based on this calculation: 
the first simply prop agates the error in iV(SII), while the second includes the e rrors in 
the s olar abundances ( Asplund et al.ll2005l ) and the natural variation of depletions (Jenkins 
2004] ). Figure [4] demonstrates that S II is a reasonable proxy for C II, where 12/20 (60%) 
agree within 2a, and those that disagree by >2a are all overestimated by measuring the 
strongly saturated C II resonance line directly. Forcing A^(CII) = A^(CIIsii) in fitting C II 
and C ir produces viable fits to the data that are indistinguishable from those shown in 
Figures [la] and [lb] The discrepancy between the iV(CII) and A^(CIIsii) originates from the 
well-known problem when on the fiat part of the curve-of-growth, that very small changes 
in the Doppler width can produce very large changes in the column density. Despite the 
complication of using the ISM average depletions for carbon and sulfur, estimating the C II 
column density is likely more accurately achieved by using S II. 



3.5. Mg I and Mg II 

Spectra of three (r/ UMa, a Gru, and G191-B2B) of the 13 sight lines with C IF 
absorption also contain the Mg I and Mg II lines. We fit the optically thin Mg I line 
separately and the two marginally saturated Mg II lines simultaneously. The spectra and 
fits are shown in Fig [3] and the fit parameters are listed in Table 15] Due to the strength of 
the Mg II lines and relative weakness of the Mg I, not all Mg II components are detected 
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in Mg I. Although those that are, have consistent velocities and Doppler widths, indicating 
both ions are likely part of the same collections of gas. 



3.6. Individual Sight Lines 

a Aur. — The line of si ght toward a Aur is well k nown to exhibit a strong single interstel- 



lar absorption compone nt (Linskv et al. 



1993 



19951 ) that has been identified with the Local 



Interstellar Cloud (LIC: iLallement et al.lll995l ). These characteristics, complemented by the 
star's proximity, contributed to it being t he first LISM sight line with a Uf. measurement based 
on absorption detected in C II and C II* (IWood &: Liiiskvlll997l ). The first fit plotted in Fig- 
ure[Ta|is for a Aur. As noted by lWood fc LinskyI (119971 ). who analyzed the same data, the low 
signal-to-noise and the saturation of the resonance line make the measurement of a precise 
Doppler width difficult. For this reason, we use the independently determined LISM temper- 
ature a nd turbulent velocity along this line of sight to force the Doppler width to be 3.48]'^Q;Jg 



km s ^ (JRedfield fc Linskyll2004bl ). Our measure d central velocity of t^ = 20.78 ±0.28 km s ^ 
agrees with the 20.2 ± 1.5 km s^^ measured by IWood &: LinskyI (119971 ). Also, we measure 
column densities of logiV(C II)= 14.67+o:22 and logA^(C II*)= 12.62 ±0.07, which agree with 
their logA^(C II)= 14.8 ± 0.3 and logA^(C 11*)= 12.64 ± 0.07, respectively. 

r\ UMa. — The rj UMa sight line shows two components. The Doppler width of the 
absorption observed toward this ta rget is constrained based on independent temperature 
and turbulent velocity information (JRedfield fc LinskyI l2004bl ). For the first component, b 



and for the secon d comp onent, b is fixed at 5.60]^o ii k'^ ^ 



-1 



is fixed at 3.76+E5:?? km s'^ 
All C II fit parameters from iRedfield fc LinskyI (l2004al ) agree with our measurements. Both 
components are detected in Mg I and Mg II, although the redward component is very weak 
and has a poorly determined central velocity and Doppler width. Indeed, the blueward 
(~-2 km s~^) component is consistent in velocity for both carbon and magnesium, but the 
measured velocity of the redwa rd componen t is si gnificantly discrepant between the two, 
although the errors are large. iFrisch et al.l (120061 ) provided an in-depth analysis of this 
particular sight line, including a measurement of ne ~ 0.1 cm~^. 

a Gru. — The spectra of a Gru were taken at the lowest resolution allowed in our survey. 
Based on this observation alone, it is difficult to identify the appropriate number of absorp- 
tion components, but on the merit of past high spectral resolution observations of other ions 
(e.g., Mg II, Fe II), three inte rstellar components have been determined along this line of 
sight (JRedfield fc LinskyI |2002| ). We fit three components accor dingly and use the same dif- 
ference in central velocities between the three components as in iRedfield fc LinskyI (120021 ) to 
constrain our C II fit. We use the wavelength solution of Fe II in order to calibrate the C II 
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spectrum, which had no wavelength cahbration image at the same wavelength. The errors are 
large for Doppler width and column density due to the blending of the three components. 
Our rn easurements are consistent, though more reliable, than those by iRedfield fc Linsky 
( l2004al ). in w hich the C II reson a nce li ne was fit alone. Our Mg II measurements agree well 
with those of iRedfield fc Linskvl J2002h . 



WD 0050-332.— Vsing low-resoluti on (R ~10,000 : Av ~ 30 km s'^) International Ul 



traviolet Explorer (lUE) observations, iHolberg et al.l (119981 ) find interstellar C II toward 
WD 0050-332 at 10.69 ± 3.20 km s~^. At higher resolution, our C II and S II observations 
resolve this feature into two components at ~6.3 and ~12.6 km s~^. In agreement with 
ou r claim of two compo nents, two components have also been observed in N I and Si II 
by lOliveira et al.l (120051 ) at ~6.8 and ~16.8 km s~^. In C 11*, however, we observe only 



the stronger ~6.3 km s ^ component. Tentative detections of highly ionize d components 



associated with a circumste llar shell have been observed toward this target (IHolberg et al. 
19981 : iBannister et al.l[2003[ ). 



EX Hya and IX V el. — Both targets have been scrutinized for observational evidence 
for circumbinary disks (JBelle et al.l 120041 ). In this search, which used the same datasets 
employed here, they identify absorption in C II, C II*, and S II, altho ugh the absorption 
was not identified with circumbinary absorption for either star. Instead, iBelle et al.l (120041 ) 
attributed the observed absorption to the LISM and fit the profile with a single Gaussian 
function, although particularly in the case of IX Vel, two components are clearly present. We 
fit both EX Hya and IX Vel with two components, motivated by clear asymmetries in the 
C II* and S II absorption. The column-weighted velocity of our components match fairly well 
with their estimates for EX Hya. However, the weaker C II* transition (1335.6627 A) appears 
to be used as the wavelength standard instead of the stronger transition (1335.7077 A), 
which overestimates the velocity of absorption by ~10 km s^^. This correction improves 
the agreement of their C 11* velocity with the other lines they measured and with our 
measurements. Their IX Vel estimates are significantly different from ours and are likely 



due to the poor approximation of a single Gaussian to the absorption profile. iLinnell et al. 



(120071 ) also note the presence of interstellar absorption in C II and S II in the spectrum of 
IX Vel. 

G191-B2B. — This is the third target which the Doppler width c an be constrained 
from i ndependent measurements of temperature and turbulent velocity (IRedfield fc Linsky 



4.1 01°-^? km 



and h 



o 40+0.21 

'J-^'J-0.26 



2004bl ). For the first and second components, h 
km s~^, respectively. The column density measured by IRedfield fc LinskyI (J2004al ) for the 
~17 km s~^ component, logA^(C II)= 15.70 ± 0.36, is consistent with our measurement, 
logA^(C II)= 15.42 ± 0.17. The ~6 km s~^ component is in significant disagreement, with 
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their esti mate at logiVfC II) = 16.12 ± 0.33 and ours at logA^(C II)= 14.8+?:2. The measure- 
ments by iRedfield fc LinskyI (J2004aJ) were obtained from fitting the saturated C II resonance 
hne only, which typically overestimates the column density. The current fits are more reliable 
since both C II and C II* profiles are used simultaneously, which more tightly constrain the fit 
parameters. We measure logiV(C II*)= 12.28± 0.09 for the second com ponent, which matches 
a velocity believed to be attributed to the LIC (JLallement et al.lll995l ) in this direction. Lying 
angularly close (~2.5°) and believed to traverse the LIC is the line of sight towards a Aur. 
We measure its column density to be slightly greater with logA^(C II*)= 12.62 ± 0.07. Since 
the distance to a Aur is 20% that to G191-B2B, it is clear that the LIC does not extend 
beyond ~13 pc in that direction. Mg I, Mg II, and S II are also measured along this line of 
sight using two components. The velocities agree between magnesium and sulfur, although 
they are slightly larger than those measured with carbon. Indeed, th e se C II and C II* fits 
differ in velocity even from the C II fits alone from iRedfield &: LinskyI (J2004al ). Again, since 
the current measurements are also constrained by the optically thin excited line, these mea- 
surements are more reliable. Why the G191-B2B absolute radial velocities disagree among 
different ions is unknown, but may be explained by systematic wavelength calibration issues, 
as the difference in radial velocity between the components agrees for all ions, and the spectra 
used for different ions a r e coad ditions of different individual observations at various grating 
settings. iLemoine et al.l (119961 ) also measure LISM absorption for two dominant components 
at ~9.9 and ~20.6 km s~^ and derive roughly similar column densities for Mg II and C II. 



WD O^g^+O^J.— iDupree fc Raymond! (119821 ) identified interstellar absorption in O I 



Si II, and N I while suggesting the presence of a Stromgren sphere due to absorption observed 
in C IV. There ha ve also been continued discuss i ons of high ionization circumstellar feature s 
for this star ( e. g.. IVennes fc ThorstensenI 1 1994 : iHolberg et al.l Il998l : iBannister et al.l |2003|) . 
Holberg et aU Jl998l ) observe LISM absorption of C I I at t; = -4 .28 ± 2.78 km s"^ with 



low resolution lUE spectra. Using the same datasets as lVennes et al.l (120001 ) . we resolve two 
components in both C II and S II at ~3.9 and ~17.0 km s~^, which match their ~3.1 and 
~17.6 km s~^ components very well. We find evidenc e of the ^17.0 km s~^ component in 
S II, most notably in the strong 1259.5 A line, whereas IVennes et al.l ( 120001 ) only fit the ~3.9 
km s~^ component. We observe, just as they do, only the ~3.9 km s~^ component in the 
C ir excited line. 



WD 2309+105.- 



Holberg et al.l (119981 ) use low-resolution lUE spectra to identify one 
LISM component in C 11* aX v = —8.27 ± 3.19 km s^^ and in all three S II lines at velocities 
ranging fror n -14.5 to -6.7 k m s~\ Two components are seen in high resolution if ST spectra 
analyzed by lOliveira et al.l (120031 ) in N I, S II, and Si II, while only the redward component 
is detected in C 11*. Using the same dataset, we concur, and derive nearly identical fit 
parameters. Our central velocities are slightly different; ours are measured at ~-8.5 ± 
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1.2 and ~1.3 ± 1.8 km s~^ and theirs at ~-9.8 and ~-0.6 km s^^, but the differences 
are minimal given the large errors and almost identical differential velocity. Only the ^-8.5 
km s~ ^ component is detected in C 11* in agreement with the /t/^i? detection by lHolberg et al. 
Jl998[ ). 



WD 121 0+533. — Although on ly one component could be discerned in the low resolution 
lUE spectra, iHolberg et al.l (119981 ) find C II absorption at f = —7.32 ± 3.21 km s~^ and 



C ir absorption at f = —12.14 ± 5.92 km s~^. We detect three components in C II, while 
only the two redward components are seen in S II, at velocities of ~-23.4, ~-9.3, and ~-3.0 
km s~^. The S II 1253.8 A line suffers from considerable contamination, although the two 
components are well characterized in the other two S II lines. 

p Lup. — The two components characterizing the excited line of p Lup are similar in 
Doppler width and in column density, lending difficulty to distinguishing the two despite 
a ~7 km s~^ separation. Fortunately, the asymmetric S II absorption profile allows us to 
discern two LISM cor nponents. We use v e lociti es determined from these measurements to 
constrain our C II fits. I Welsh fc LallementI ( 20051 ) observe two components in Na I (5890.0 A) 
and three components in Fe II (1608.5 A), S II (1253,1259 A), and Al II (1670.8 A). Ahhough 
we only find evidence for two components in C II and S II, our compon ents agree fairly well, 
with velocities ~-16.1 and ~-9.1 km s~^. I Welsh fc LallementI (120051 ) present a thorough 
discussion of cloud distances and identify the ~-9.1 km s~^ absorption with material at ~90 
pc while identifying the ~-16.1 km s~^ absorption with material at the neutral boundary of 
the Local Bubble in the direction toward p Lup. 



4. Electron Density 



4.1. Estimation Based On C II* 



Table [6] lists the measured electron densities, Ue, along the lines of sight towards the 
13 targets analyzed. The values hsted in column four of this table are calculated using a 
method which compares the column densities of the res onance and excited lines of C II. 
Our use of thi s meth od is similar to that imple mented bylSpitzer fc Fitzpatrickl (Il993l ) and 
Oliveira et al.l (120031 ) but most parallels that of IWood fc LinskyI (119971 ). 



The C II resonance absorption line at 1334.5323 A corresponds to the transition from 
the ground state (J = 1/2), while the C II* excited absorption lines at 1335.6627 A and 
1335.7077 A correspond to the transition from the excited state of the fine-structure doublet 
(J = 3/2). Collisions with electrons are responsible for populating the excited state, and 
hence the ratio of the column densities of the two lines is proportional to the electron density. 
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For a detail ed discussion of how the fine structure of absorption fines can be used to determine 
density, see lSaficafi fc Wolj JlQGsl ). 



Tfie fofiowing relation is derived from tfiermal equilibrium between collisional excitation 
of tfie J = 3/2 state and radiative de-excitation: 



iV(CIP) _ n,Ci2{T) 
N{CU) A21 



(3) 



Tfie effect of collisional de-excitation at tfiese densities and temperatures is negligible and 
tfierefore not included in tfie equation above. N[C II) and A^(C II*) are tfie column densities 
of tfie resonance and excited lines, respectively. Tfie calculation of tfie electron density 
using S II as a proxy for C II, simply replaces A^(C II) witfi A^(C IIsn), as derived from 
Equation [2l Tfie electron densities based on A^(C IIsn) are also listed in Table [6] (column 
five). Tfie radiative de-exci tation rate coefficient is A21 = 2.29 x 10~^ s~^, as listed by 
Nussbaumer fc Storeyl (1198 ll ). Tfie collision rate coefficient can be expressed is cgs units as 



Ci2iT) 



8.63 X lO-^fi 



12 



9iT'' 



exp 



E 



12 



kT 



(4) 



wfiere tfie statistical weigfit of tfie ground state gi = 2, and tfie energy of tfie transition 
E12 = 1.31 X 10~^^ ergs. Since tfie collision s trengtfi Q^9, fias a very wea k temperature 
dependence, for all targets we let fli2 = 2.81 (JHayes fc Nussbaumerl Il984l). calculated at 
a temp erature of 7000 K, very similar to tfie LISM average of 6680 K fJRedfield &: Linsky 
2004bh . 



4.2. Estimation Based On Mg I and Mg II 

Tfiree of our targets fiave botfi Mg I and Mg II LISM absorption measur ements (see Ta- 
ble [5l) . We estimat e tfie e lectron density by following tfie procedure detailed by lLallement fc Ferlet 
( 119971 ) and iFriscfil (119941 ). Assuming equilibrium, tfie electron density can be estimated by 
formulating tfie balance of ionization and recombination between neutral and singly ionized 
magnesium: 

A^(Mgll) ^ r + ne^ex , . 

A^(Mgl) riea ' ^ ' 

wfiere F is tfie pfiotoionization rate of Mg I, cXex is tfie formation rate of Mg II based on cfiarge 
excfiange, a is tfie total recombination rate, and we fiave assumed tfiat A^(Mg II)/A^(Mg I) = 



n(Mg II)/n(Mg I) and n(H II) 



He 



We u se a pfiotoionization rate of F 



6.1 X 10^1 s-^ 



(jSofia &: Jenkins! Il998l : I Jenkins et al.ll2000al ). Altfiougfi pfiotoionization typically dominates 
over cfiarge excfiange, at fiigfi enougfi densities, tfie cfiarge excfiange contribution can be 
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significant. We have used the relationship derived by I Allan et aLl (119881). ( Xex = 1.74 x 



10 ^exp(— 22100/T) cm^ s ^. The total recom bination rate, similar to 



sum of the radiative recomb ination rate from lAldrqvandi fc Peauignot 



ll973h and d ielec- 
tric recombination rate from iNussbaumer fc Storey! (Il986l ) and iMazzotta et al.l (119981 ) . All 



rischl (Il994f) . is the 



recombination rates are strongly dependent on temperature. Table [6] lists the calculated 
electron densities based on the ratio of Mg I and Mg II (column six). 



4.3. Previous Measurements 



Electron density estimates have been made by other researchers for a few of the sight 

+0.060 
-0.059 



lines. Toward a Aur, our electro n density {rip = C).140l nn^Q cm ^) agrees well with the 
measure d n^ = 0.11^'^nnK cm~^ by IWood fc LinskyI (119971 ). An interstellar temperature of 



6700 K (IRedfield &: Linskvl l2004bl ) is included in our determination of the electron density. 
Vennes et al.l (120001 ) obtained an electron density estimate toward WD 0232+035 by assum- 
ing a wide range of possible b values which resulted in a density of rig = 0.36 cm~'^ for the 
interstellar component near 3.8 km s~^. This is slightly high er, but within 2(7 of our esti- 

3 and ne(CIIsii) = 0.21^°;^^ cm'^. lohveira et all J2003h provide 



h0.17 
-0.6 



cm 



mates, He = 0.10'^^ 

an estimate of n^ = 0.1 ± 0.01 cm~^ for the -9.7 km s~^ component toward WD 2309+105, 
based on using S II as a proxy of C II and assuming a temperature of 8000 K. Our measure- 
ment using the same technique, assuming the LISM average temperature of 6680 K, matches 
quite well, 



n. 



0.084+°:Si cm-3. 



Discussion 



Figures [5H7| summarize the results of our electron density measurements. Figure [5] shows 
the measured electron density as a function of the C 11* column density. Electron density 
estimates using both the saturated C II resonance line and S II as a proxy for C II are shown. 
The measurements using S II are more precise, even when we include systematic errors of 
the variation of depletion, than the C II measurements since S II is optically thin. Figure [6] 
displays the same distribution of data points, except instead of electron density, we plot 
the resonance line column density versus the excited line column density. Overplotted are 
constant-density contours for the mean temperature of the LISM. Both plots show that the 
LISM electron density measurements are relatively tightly log-normal distributed about the 
unweighted mean value of ~0.1 cm"^. Figure [7] further emphasizes the tight distribution 
of electron density measurements about the mean. In particular, by using optically thin 
S II lines as a proxy for C II, we avoid systematically high column density measurements 
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commonly derived from saturated line profiles, which ultimately lead to very low electron 
density determinations. For this reason, in the electron density distribution derived from 
S II, we lose the low electron density tail evident with the C II measurements, and the 
LISM average tightens to an unweighted mean of O.ll^QQg cm~^, where the la errors are 
the dispersion about the mean value. All electron densities derived using S II as a proxy 
range from 0.07-0.80 cm~^. The distribution of ne(C IIsii) measurements matches well the 
ne(C II) distribution, but it is significantly tighter, since the gross overestimates of the C II 
column density are avoided. This is a strong endorsement for using S II as a proxy for C II, 
in order to make more accurate electron density measurements. 

Figure [8] shows the solutions of electron density as a function of temperature for the 
C ir/C II technique, which does not require ionization equilibrium and has a very modest 
temperature dependence, and for the Mg I/Mg II technique, which does require ionization 
equilibr ium and has a v ery st rong temperature dependence. This type of comparison was 



used by iGry fc Jenkins! (1200 ll ) to measure the properties along the sight line toward e CMa. 
Also shown in Figure [8], for several of the sight lines, is an independent measurement of 
the temperature of the gas based on the widths of absorption lines of different atomic mass 



(JRedfield fc LinskyI l2004bl ). Ideally, all three measurements should converge at the same 



temperature and provide a single estimate of the electron density. 



5.1. LISM Cloud Properties 

The last column of Table [6] indicates the interstellar cloud(s) by which the observed 
components may be identified. The criteria for identification include agree ment with the pre- 



dicte d projected velocity and spatial distribution of a given LISM cloud (JRedfield fc Linsky 



20081 ). Components for some of our targets (a Aur, a Gru, and G191-B2B for the LIC; 



r/ UMa for the NGP cloud; G191- B2B for the Hyades clou d) had been previously identified 



and assigned to specific clouds by lRedfield fc LinskyI (120081 ). Many of the remaining compo- 



nents could also be identified with known LISM clouds (i.e., the observed radial velocity of 
absorption is within 3o" of the predicted projected velocity, and the line of sight traverses the 
spatial distribution of the cloud). Those clouds listed in parentheses also agree in velocity, 
but the sight line only passes near (within 20 degrees of) the given cloud boundary. Given 
the similarities of the LISM cloud velocities, it is not always easy to uniquely identify cloud 
membership, although of the 23 C 11* components, 11 can be firmly identified with a specific 
LISM cloud. 

The Local Interstellar Cloud (LIC) is the collection of gas that dominates LISM ab- 
sorption line observations because it is detected in a large fraction of the sky. The Sun is 



-17- 



currently located just outside of the LIC (JRedfield fc Linskvll2000l . l2008l ). While LIC mate- 
rial has likely surrounded the Sun for the last ~100,000 years, given the relative motion and 
location of the LIC and the Sun, the solar system has moved, or will very shortly move, into 
a different interstellar environment. However, the proximity of the LIC means that many 
LISM sight lines will probe this material. Indeed, 7 of the 11 identified cloud members are of 
the LIC. The LIC Ue measurements are shown by the shaded histogram in Figure [71 where all 
sight lines give a similar measurement. The ne(C IIsii) measurements are used to calculate 
the mean when available, while ne(C II) is used for a Aur and ne(Mg II/Mg I) is used for 
a Gru. The weighted mean is rae(LIC) = 0.12 ± 0.04 cm~^. Dramatic variation in electron 
density within the LIC is not observed. However, more sight lines would be required to do 
a more involved investigation of intracloud variability. 

Other electron density components are assigned to various other clouds, but none have 
as many measurements as the LIC. The North Galactic Pole (NGP) Cloud has two He 
measurements which agree quite will with each other at ~0.09 cm~^ and are not too different 
than the LIC measurement. Likewise, the Gem Cloud is observed in the second component 
of EX Hya (rig ~ 0.15) and is similar to the LIC measurement. This is the first strong 
evidence that LISM clouds that are dynamically distinct from the LIC nonetheless have 
a similar electron density. However, not all LISM clouds share similar electron density 
properties. The Hyades Cloud measurement observed in the first component of G191-B2B is 
consistently measured at a high electron density (~0.5) in all three techniques. The Hyades 
Cloud appears t o be a decelerated cloud at the leading edge of the platoon of LISM clouds 
beyond the LIC jRedfield fc LinskvlboOliboOsh . It is therefore closer to G191-B2B than the 
LIC and without any obvious LISM clouds between it and G191-B2B to shield the ionizing 
radiation. The Hyades Cloud appears to have an enhanced electron density due to the 
photoionization of G191-B2B, while also shielding other LISM clouds, such as the LIC, from 
the strong ionizing radiation of G191-B2B. No obvious measurements of the second largest 
LISM cloud and likely future interstellar environment of the Sun, the G Cloud, exist in this 
sample. 



5.2. Sources of Ionization and Self-Shielding 



Almost all of our background stars are white dwarfs or early-type stars and therefore 
strong UV pho ton sources and significant contributors to the local radiation field. Indeed, 
Vallergal (119981 ) provided an inventory of 54 strong extreme-UV (EUV) stars that largely 
determine the local radiation field. The local EUV field is dominated by the B star e CMa 
and by three white dwarfs (Feige 24, HZ 43, and G191-B2B). Both Feige 24 and G191- 
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B2B are identified here as having detectable C II*, while e CMa is a relatively distant star 
(~130 pc) and not included in our search sample, and HZ 43 did not have any detectable 
C II*, although this may be due to the significantly lower amount of LISM material in that 
direction. Three more of our 13 sight lines (a Aur, Wp0050- 332, and WD2309+105) are also 



included in the radiation field inventory by IVallergal (Il998l ). and we expect that the other 



white dwarfs and early-type stars that make up our remaining sight lines are also significant 
contributors to the local ionizing radiation field. 

In the left plot of Figure M, the electron density measur ements are sho wn as a function 



of the angle from e CMa. This is a test that was suggested in IVallergal (119981 ). Unfortunately, 
most sight lines cluster around 90 degrees from the direction of e CMa, but nonetheless, no 
clear correlation is detected. This constrasts with a possible ionization gradient toward e 
CMa, based on hydrogen and helium c olumn densities rn easured in the extreme ultraviolet 



toward many nearby white dwarf stars (IWolff et al.lll999l ). Our measurements indicate that 



e CMa does not singularly dominate the ionization structure of the LISM, and other con- 
tributors, such as the other nearby WDs and early type stars, may be significant sources of 
the radiation field that dictates the electron density in the LISM. 

The right plot of Figure |9] demonstrates a possible test of self-shielding in LISM clouds. 
The figure plots the difference of electron densities along the same sight line as a function 
of the total column density along the line of sight. Those with low columns should show 
evidence of less shielding, while large columns will provide significant shielding and a more 
dramatic difference in electron density measurements. No such correlation is seen, indicating 
that for the entire LISM sample, there does not appear to be pervasive shielding of ionizing 
radiation along each specific line of sight. This diagnostic assumes that for each sight line, 
the dominant ionization source is the background star itself, which will not necessarily be the 
case. A more sophisticated three-dimensional morphological model of the LISM that includes 
the effects of known ionizing sources is needed to accurately determine if the observed electron 
density is consistent with shielding of LISM clouds. However, isolated sight lines of strong 
ionization sources may be used to investigate the impact of shielding on electron density 
measurements. 

As mentioned in the previous section, the G191-B2B sight line provides a possible exam- 
ple of self-shielding. Because many of our background targets are significant EUV sources, 
clouds closer to the source (i.e., more distant from the Sun) should shield the clouds farther 
from the source (i.e., nearer to the Sun). The LIC components offer excellent examples of 
this, since we know the LIC component is tracing gas closest to the Sun. The G191-B2B 
sight line, and others including the a Gru, WD0232+035, and WD1210-I-533 sight lines, 
are cases where the LIC component electron density is less than the other observed compo- 
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nent, presumably probing material from a more dista nt cloud that lie s close r to the ionizing 
source (assuming it is the background star observed). iGry &: Jenkins! (1200 ll ) also noted this 
phenomena in looking at the electron densities toward the most dominant ionization source, 
e CMa. They identified two absorbers with local clouds, the LIC and Blue Cloud, which had 
low electron densities, He =0.08-0.17 and 0.016-0.088 cm~'^, respectively, while the distant 
absorber, which was located closer to e CMa, had a significantly higher electron density. 
He =0.18-0.28 cm~^. Several counter-examples, such as GD659, and 74 Psc A and B, are 
also observed, where the LIC does not have a significantly lower electron density. However, 
these are among the weakest ionizing sources in our sample. The 74 Psc stars are inactive 
and relatively cool A stars, while GD659 is t he coolest (~35000 K) of the singleton white 



dwar fs in our sample, which are all >50000 K fJLajoie &: Bergeronll2007l : iHolberg &: Bergeron 
2006h . 



5.3. LISM Pressure Measurements 



The electron density is also important for its implication for pressure measurements 
of LISM material. This is a particularly critical issue due to the apparent disparity of 
the warm LISM clouds, such as those observed by C II* abso rption, which have pressures. 



P/k = riT ^ 3000 K cm-^ (e.g., iRedfieldl l2006l : I Jenkins! l2002[ ) . and the hot, tenuous Local 



Bubble gas that sur rounds the warm clouds, which has pressures ~10000 K cm ^ (e.g. 



Snowden et al.! !l990! ) . Recent evidence of soft X-ray emission at the heliosphere due to 
charge exchange between the solar wind and incoming LISM appears to contribute to the 
soft X-ray emiss ion that was previously fully attributed to the hot Local Bubble gas (e.g., 
Lallement! !2004! ). A revised inventory of soft X-ray emission may lower the temperature 
and/or density of the hot gas and reduce the pressure discrepancy. 

The top axis of Figure [71 shows the range of measured pressures that we obtain from 
our electron density measurements. This calculation assumes temperature is constant for all 
sight lines, and a simple photoionization relationship between the electron density (rig) and 
the neutral hydrogen density (nni). In the case of temperature, we know this assumption is 
not completely valid sinc e we see some variation about the mean LISM value (T = 6680 K 
Redfield fc Linsky!!2004a! !b!). However, since independent temperature measurements are not 
available for the majority of sight lines studied here, and since the dispersion about the mean 
temperature is not high, it is a reasonabl e initial assumption. We assume riui = ''^e'^(H)/r(H) 
(see Equation 7 in !Sofia fc Jenkins!!l998! ). The balance of the recombination rate (a) and the 
ionization rate (F) are assumed to be constant in the LISM and we calibrate this quantity 
such that the LISM average electron density (rig) is consistent with the LISM measurement 
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of nH T = 0.222 cm ^, based on in situ measurements of riuci = 0.0151 cm ^ ( IGloeckler et al. 
2004). and the H I to He I column density ratio observed toward nearby white dwarfs 



( JDupuis et al.l 119951 ). We assume Up = n^ and riHe = O.lnn. This calculation breaks down 
at very low densities (rig < 0.007 cm~^), where the derived hydrogen density is less than 
the minimum value allowed from observations of nearby stars, measured by dividing the ob- 
served hydrogen column d ensity by the distance to the background source (JRedfield fc Linsky 
2004al : iLinsky et al.ll2000l ). At high densities (ng > 0.4 cm~^), our assumption of a constant 
ionization rate fails severely because the resulting column density thro ugh a typical LISM 
cloud (A^Hi > 10^^'^) leads to significant shielding of ionizing ratiation (jjenkinsll2004l ). The 
top axis of Figure [7] is only printed for the range of densities for which our calculation is 
reasonable. The distribution of electron densities then translates into an unweighted mean 
P/k = 3300l^gQQ K cm~^, consistent wi th the r a nge o f values determined for other nearby 
stars using excited transitions of C I by I Jenkins! (120021 ). 



Conclusions 



We analyze high spectral resolution observations of LISM absorption in order to sur- 
vey the electron density in nearby interstellar material. These measurements should pro- 
vide important constraints on the ionization and abundance p atterns of the LISM (e.g., 
Slavin &: Frischll2002l : ISofia &: Jenkinslll998l: I Jenkins et al.ll2000bf). as w ell as on evolutionary 



models of all phases of the LISM ( Breitschwerdt fc de Avillea l2006l ) . A summary of 
results is as follows: 



our 



1. We searched the entire HST spectroscopic database of nearby stars (<100 pc) for 
detections of C 11*. Of the ~417 total nearby sight lines, we find 13 that show C 11* 
absorption in 23 different velocity components. The vast majority of these detections 
are new. 

2. Using the C 11* to C II ratio, we infer the electron density. To increase the accuracy of 
our results, particularly in terms of measuring the column density of the saturated C II 
resonance line, we employ three analytical strategies: (a) simultaneously fitting both 
the C ir and C II profiles, allowing the optical thin C II* line to constrain the line 
width, (b) using independently derived temperatures from comparison of line widths 
to constrain the acceptable range of line widths for C II, and (c) using easily measured 
S II column densities as a proxy for C II column density. 

3. The distribution of electron densities based on using S II as a proxy for C II is similar 
to the distribution based on carbon alone, while significantly tighter. This is a promis- 
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ing technique to avoid grossly overestimating the C II column density based on the 
saturated line profile. 

4. We find the distribution of measured LISM electron densities (ue) is consistent with a 
log-normal profile, with a mean (unweighted) value of ne(C IIsn) = 0.111o;q5 cm^^. 

5. We assign individual velocity components to specific LISM clouds based on kinematical 
and spatial properties. In particular, the LIC is probed by seven different sight lines, 
which all give roughly identical electron density measurements. The weighted mean 
value for the LIC is Ug = 0.12 ± 0.04 cm'^ 

6. Two clouds, the NGP and Gem clouds, show similar electron density properties as 
the LIC. The Hyades Cloud, a decelerated cloud at the leading edge of the platoon of 
LISM clouds, has a significantly higher electron density than the LIC. Observed toward 
G191-B2B, the high electron density may be caused by the lack of shielding from such 
a strong radiation source. 

7. Almost all of our background sources are significant ionizing sources that may infiuence 
the ionization structure and thereby the electron density of the gas along the line of 
sight. We do not find evidence that the ionization structure of the LISM is dominated 
by a single source, namely e CMa. 

8. We see evidence of more distant clouds (i.e., those closest to the ionizing sources) 
shielding nearer clouds (i.e., those farthest from the ionizing sources). In several ex- 
amples, the LIC component which is known to be farthest from the ionizing source 
has a lower electron density than the component along the same line of sight that is 
nearer to the radiation source. Although counter-examples exist, they are toward the 
weakest radiation sources in our sample. 

9. The range in electron density is used to estimate the range of pressures that may 
be found in warm LISM clouds. Given simple assumptions, the measured electron 
densities correspond to an unweighted mean pressure P/k = 3300ti9oo -^ cm~^. 
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Table 1. Stellar Parameters 



HD 


Other 


Other 


Spectral 


/ 


b 


Distance'^ 


# 


Name 


Name 


Type 


(deg) 


(deg) 


(pc) 


34029 


a Aur 


Capella 


G0III+G8III 


162.58 


+4.566 


12.94 ± 0.15 


120315 


7] UMa 


Alcaid 


B3V 


100.69 


+65.32 


30.87 ± 0.71 


209952 


a Gru 


Alnair 


B7IV 


349.99 


-52.47 


31.10 ±0.79 




WD 0050-332 


GD659 


DAl^ 


299.14 


-84.11 


-58 




EXHya 




M5/M6+DA 


303.18 


+33.62 


64.5 ± 1.2 


6457 


74 Psc B 




AOV 


127.34 


-41.28 


68.2 ± 8.0 




WD 0501+527 


G191-B2B 


DAl^ 


155.95 


+7.099 


69 ±15 


6456 


74 Psc A 




AlV 


127.34 


-41.27 


73.2 ± 7.9 




WD 0232+035 


Feige 24 


DAZQOl"^ 


165.96 


-50.26 


74 ± 20 




WD 2309+105 


GD246 


DAl"^ 


87.262 


-45.11 


-79 




WD 1210+533 




DAO'^ 


135.61 


+63.11 


-87 


128345 


p Lup 




B5V 


320.13 


+9.857 


95.1 ± 6.4 




IX Vel 




B8V+DA 


264.92 


-7.890 


96.3 ± 9.1 



Note. — All values from SIMBAD unless otherwise noted. 



'^AU distances are Hipparcos distances flPerrvman et al.lll99'jt) exc ept for WD 0050-332 
and WD 2309+105 w hose distances a re fro m Ivennes et all (|l997l ) and WD 1210+533 
whose distance is from iHolberg et al.l (119981 ). 



^White dwarf spectral types taken from lMcCook fc Sioru (119991 ). 
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Table 2. Observational Parameters 











Spectral 




Exposure 




Other 


Other 






Range 


Resolution 


Time" 




Name 


Name 


Instrument 


Grating 


(A) 


(A/AA) 


M 


Dataset 


CK Aur 


Capclla 


GHRS 


ECH-A 


1331-1338 


100000 


761.6 


Z2UW030BT 


r, UMa 


Alcaid 


GHRS 


ECH-A 


1332-1339 


100000 


1171.4 


Z3CL020ET 






GHRS 


ECH-B 


2790-2805 


100000 


108. 8 


Z3CL030HT 






GHRS 


ECH-B 


2845-2859 


100000 


652.8 


Z3CL030ST 
Z3CL0307T 


a Gru 


Alnair 


GHRS 


G160M 


1309-1345 


20000 


230.4 


Z1720109T 






GHRS 


ECH-B 


2791-2806 


100000 


108. 8 


Z1720209T 






GHRS 


ECH-B 


2847-2861 


100000 


217.6 


Z172020AT 


WD 0050-332 


GD 659 


STIS 


E140H 


1170-1372 


114000 


4134 


O4G101010 
O4G101020 


EX Hya 




STIS 


E140M 


1140-1735 


45800 


15200 


068301010 
068301020 
068301030 
068302010 
068302020 
068302030 


74 Psc B 


HR 311 


STIS 


E140H 


1170-1372 


114000 


2768 


O56L02010 
O56L02020 
O56L02030 


WD 0501+527 


G191-B2B 


STIS 


E140H 


1170-1517*' 


114000 


5623 


O57U01020 
O6HB10040 
O6HB10050 
O6HB10060 
O6HBl()07() 
O6HBl()080 
O6HB10090 






STIS 


E230H 


2624-3095*' 


114000 


6013 


O6HB30080 
O6HB30090 
O6HB300B0 
O6HB300C0 
O6HB300D0 
O6HB300E0 


74 Psc A 


HR, 310 


STIS 


E140H 


1170-1372 


114000 


4128 


O56L01010 
O56L01020 
O56Lai030 
O56L51020 


WD 0232+035 


Feigc 24 


STIS 


E140M 


1150-1735 


45800 


4176 


O4G701010 
O4G702010 


WD 2309+105 


GD 246 


STIS 


E140H 


1170-1372 


114000 


2420 


O4G102020 


WD 1210+533 




STIS 


E140M 


1140-1735 


45800 


8371 


O5F203010 
O5F203020 
O5F204010 


p Lup 


HR 5453 


STIS 


E140M 


1150-1735 


45800 


1900 


O8S602010 


IX Vol 




STIS 


E140M 


1140-1735 


45800 


5250 


O5BI01010 
056102010 
O5BI03010 



^For targets with multiple datasets, the exposure time listed is the sum of the exposure times for each dataset. 
Spectral range encompassed by multiple datasets. 
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Table 3. Fit Parameters for C II ISM Velocity Components 



Other 


Other 


Comp. 


V 


b 


logAf(C 11)'' 


logAf(Ciisii) 


logAf(Cn*)t= 


Name 


Name 


# 


(km s-i) 


(km s-l) 


log (cm-2) 


log (cm-2) 


log (cm-2) 


a Aur 


Capella 


1 


20.78 ± 0.28 


n 40+0. 15c 


14.67l»:i^ 




12.62«:«^ 


r? UMa 


Alcaid 


1 


-2.09 ± 0.24 


3.76±H?^ 


14.43±°:- 




12.39«:«^ 






2 


9.72 ± 0.95 


5.60l«:??^ 


-,0 -19+0.06 




ll-60«:g^ 


a Gru 


Alnair 


1 


-22.50 ± o.eo'* 


3.35 ± 0.73 


14.98±°-^ 




12.69«;}* 






2 


-12.63 ± 0.60<* 


3.49 ± 0.51 


14.52lSi« 




12.73tO;l? 






3 


-8.32 ± 0.60^ 


3.9 ± 2.2 


umtl%t 




12 17+0-48 
J^^-t'-0.48 


WD 0050-332 


GD 659 


1 


6.17 ± 0.38 


4.31 ± 0.77 


.rr ,,+0.32 


1/1 sns+0-065 +0.216 
i4.»US_Q Qgg _o.216 


-19 v-l +0.05 
t^-'t_o.05 






2 


12.3 ± 1.0 


4.67 ± 0.67 


15.12+«-29 


14 nfi9+0-240 +0.317 
i4.UD^_Q 35Q _o.406 


<11.9 


EXHya 




1 


-15.2 ± 1.2 


5.03 ± 0.72 


16.64t°;f, 


ic 909+0.120 +0.239 
ia.Z.3Z_g J20 -0.239 


i3.oit°:lg 






2 


-6.2 ± 2.3 


6.1 ± 1.9 


14.75l«:^^ 


14 7Q9+0.290 +0.356 
i*. /»z_g 290 -0.356 


i2.78t°:r. 


74 Psc B 


HR311 


1 


-5.75 ± 0.50 


4.00 ± 0.93 


i&m+J,\l 


It; A4A+0-017 +0.207 


10 07+0.02 

t3.UC_o.o2 






2 


10.47 ± 0.70 


4.26 ± 0.44 


l^A^tlf, 


14 744+0.052 +0.213 

14. (44_g gj2 -0.213 


19 09+003 


WD 0501+527 


G191-B2B 


1 


5.98 ± 0.17 


4 ,f,+0.29c 


14.8+;i 


14 q4c+0.087 +0.224 
i4.ci40_Q ogy -0.224 


13 14+0.02 
t3.14_o 02 






2 


16.98 ± 0.60 


-0+0. 21c 


15-421^:^7 


14 Kc:n+0-052 +0.213 
i4.00U_o 052 -0.213 


12.28«;g^ 


74 Psc A 


HR 310 


1 


-5.97 ± 0.52 


3.82 ± 0.59 


16.42±«-1 


ie nQQ + 0.038 +0.210 

ia.ciyy_o 041 -0.211 


13.07«;°| 






2 


10.29 ± 0.76 


4.75 ± 0.73 


i5.04i[;:| 


14 7C9+0.110 +0.234 
i4./'b2_() jjo -0.234 


12.96«;°1 


WD 0232+035 


Feige 24 


1 


3.81 ± 0.47 


4.36 ± 0.53 


15.18_q;42 


14 O77+0.072 +0.219 
14.0 r (-0.087 -0.224 


i3.oo«:«^ 






2 


17.4 ± 4.6 


3.1 ± 1.4 


15-96lS:S 


1-? Qi 9+0.380 +0.432 
iJ.51Z (j 3g(, _o.432 


<11.8 


WD 2309+105 


GD 246 


1 


-9.7 ± 2.5° 


4.20 ± 0.34 


ic 1C+0.30 
itl.ib_(j3o 


It: 091+0.035 +0.209 
iO.OZi_g 035 _o.209 


13.05l°;°3 






2 


0.3 ± 2.5^= 


5.10 ± 0.74 


14.18t0-g 


14 71K+0.090 +0.225 
14. (lD_o.099 -0.229 


<12.3 


WD 1210+533 




1 


-23.4 ± 1.2 


3.6 ± 1.1 


13.40l«i^ 


<14.3 


<12.0 






2 


-9.0 ± 1.3 


4.3 ± 1.1 


15-50tr4^ 


1/f CQ9+0.340 +0.398 
14.D^iZ_g 250 -0.324 


13.00l°;i« 






3 


-2.3 ± 2.0 


4.8 ± 1.2 


15.44l«:«^ 


1 c; 979+0.120 +0.239 
iO.Z(Z_g J20 -0.239 


13.06t°j° 


p Lup 


HR 5453 


1 


-16.11* 


5.17 ± 0.74 


15.80t«:g 


ic; fi49+0-100 +0.229 
ia.04Z_g jQp -0.229 


13.27+°;}^ 






2 


-9.lf 


4.3 ± 1.1 


It: 1 4+0.61 
J^0-J^"t_0.61 


ic ni 9+0.310 +0.372 
iO.Ulz_g 3JQ _o.372 


13.08t°Ji 


IX Vel 




1 


4.9 ± 1.2 


5.73 ± 0.61 


14.39t«:- 


14 7nn+0-140 +0.249 
14. CUU^Q j4(, _o.249 


19 R4+006 

12.64_(,Qg 






2 


19.08 ± 0.24 


4.17 ± 0.38 


i^mtVA 


1"^ 4n9+0-024 +0.208 
10.4UZ_Q Q23 -0.208 


10 40+0.03 



'^Resonance line parameters derived from simultaneous fits. 

''Excited line parameters are weighted means of simultaneous and excited-only fits. 

'^Fixed h values based on independent temperature and turbulent velocity measurements JRedfield &: Linskvll2004lj | . 
''Fixed velocity difference between components based on measurements of interstellar Feu and Mgll IjRedfield fc Linskvll2002h . 
''Fixed velocity difference between components based on measurements of interstellar S II. 
Fixed velocities based on measurements of interstellar S II. 
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Table 4. Fit Parameters for S II ISM Velocity Components 





Other 
Name 


Comp. 

# 




Resonance Line 


a 


Other 
Name 


V 

(km s-i) 


b 

(km s-i) 


logAf 
log (cm"2) 


WD 0050-332 


GD 659 


1 


7.00 ± 0.90 


3.07 ± 0.67 


13.726 ± 0.065 






2 


15.0 ± 2.9 


2.8 ± 2.5 


12.98+0:^3^ 


EXHya 




1 


-12.7 + 2.1 


7.4 ± 2.9 


14.15 ± 0.12 






2 


-4.1 ± 4.0 


2.6 ± 1.9 


13.71 ± 0.29 


74 Psc B 


HR311 


1 


-5.59 ± 0.93 


3.48 ± 0.18 


14.362 ± 0.017 






2 


11.8 ± 1.4 


3.91 ± 0.65 


13.662 ± 0.052 


WD 0501+527 


G191-B2B 


1 


8.4 ± 1.2 


2.21 ± 0.71 


13.263 ± 0.087 






2 


20.0 ± 1.6 


4.46 ± 0.92 


13.468 ± 0.052 


74 Psc A 


HR310 


1 

2 


-6.21 ± 0.92 
12.3 ± 2.5 


4.19 ± 0.31 
4.6 ± 1.2 


,. 0-17+0.038 

13.68 ± 0.11 


WD 0232+035 


Fcigc 24 


1 


4.11 ± 0.64 


4.18 ± 0.95 


13.795tO:«r, 






2 


16.8 ± 3.8 


2.6 ± 2.3 


12.73 ± 0.38 


WD 2309+105 


GD 246 


1 


-8.3 ± 1.0 


2.42 ± 0.21 


14.239 ± 0.035 






2 


1.6 ± 1.5 


4.9 ± 2.2 


^^^ Rq4+o.o90 


WD 1210+533 




1 






<13.2 






2 


-12.0 ± 3.6 


5.0 ± 2.6 


13.55_q'25 






3 


-3.10 ± 0.94 


3.0 ± 1.0 


14.19 ± 0.12 


p Lup 


HR 5453 


1 


-16.11 ± 0.79 


5.22 ± 0.66 


14.56 ± 0.10 






2 


-9.1 ± 2.4 


4.5 ± 2.0 


13.93 ± 0.31 


IX Vel 




1 


5.3 ± 1.2 


3.8 ± 1.2 


13.618 ± 0.14 






2 


20.37 ± 0.53 


4.64 ± 0.41 


14.320t0:0^* 



^Resonance line parameters derived from simultaneous fits. 



Table 5. Fit Parameters for Mg II and Mg I ISM Velocity Components 







Other 


Comp. 




Mg II Line'^ 






Mg I Line 




Other 


V 


b 


logAT 


V 


b 


logAT 


Name 


Name 


# 


(km s^l) 


(km s^l) 


log (cm^^) 


(km s^l) 


(km s~l) 


log (cm^^) 


r, UMa 


Alcaid 


1 


-2.271 ± 0.077 


2.59 ± 0.18 


12.653 ± 0.014 


-2.36 ± 0.30 


2.59 ± 0.42 


10.153 ± 0.04 






2 


3.5 ± 1.2 


8.4 ± 3.5 


11 27+"-" 


3.5 ± 1.8 


4.8 ± 1.8 


q 27+0.13 
''■^'-0.18 


a Gru 


Alnair 


1 


-24.05 ± 0.37 


3.54 ± 0.31 


11 659+""*** 






<9.1 






2 


-14.98 ± 0.44 


3.37 ± 0.24 


13-463iS:27 


-15.09 ± 0.75 


3.43 ± 0.77 


10.99lS-^, 






3 


-7.59 ± 0.11 


1.96 ± 0.37 


12.83_o3i 


-7.6 ± 3.4 


1 Q+2-3 


io-i6-S:9o 


WD 0501+527 


G191-B2B 


1 


8.22 ± 0.23 


3.32 ± 0.34 


13 56+"-^^ 
iJ-OO_0.12 


7.94 ± 0.56 


3.38 ± 0.40 


11.193 ± 0.029 






2 


18.81 ± 0.24 


2.65 ± 0.11 


12 7-14+0.046 

li.714_„ J,g2 






<10.1 



^Resonance line parameters derived from simultaneous fits. 
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Table 6. Electron Densities 
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Comp. 
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ne(C 11311)=^ 


n,(MgII/MgI)b 


Cloud<^ 


Name 


Name 


# 


(cm-3) 


(cm-3) 


(cm-3) 




a Aur 


Capella 


1 


0.140«;g^»=* 






Lie 


7?UMa 


Alcaid 


1 

2 


0.165tnt7 
U.lt>4_Q Q27 




0.089l[5;J^^ 
0.085_o;o85 


NGP 


a Gru 


Alnair 


1 


o.08il°;gi 




<0.21 








2 


o.25i°;l« 




n 90+0.51 


(Mic, Vel) 






3 


0.5tJi 




O.wt'oil'' 


Lie 


WD 0050-332 


GD 659 


1 


0.060tO;°l 


p, ,,+0.03 +0.08d 
'^■-'-^-0.02 -0.05 






Lie, (Get) 






2 


<0.0094 


<0.11 






(Vel) 


EXHya 




1 


0.004«;SJ»^ 


n nq4+0-04 +0-074 






NGP, (Leo, G) 






2 


0.17«:f3 


n ic+o.ie +0.21 

U-iO_o.09 -0.10 






Gem, (Leo, Aur) 


74 Psc B 


HR 311 


1 


o.oi7lO;°°i 


n nfi7+0-004+o.04i 

U.UDr_Q gg4_o 026 












2 


0.036t°;°2B 


r, 1Q+0.03 +0.12d 
U-iy_o.02 -0.07 






LIG, (Hyades, Eri) 


WD 0501+527 


G191-B2B 


1 


0.30«i| 


n sn+0-18 +0-54 

u.ou_o x5 _o.33 


0-48±°:™ 


Hyades 






2 


PI f,„ +0.006 


n nsi +0-021 +o.055d 

U.U»l_o 017-0.035 


<0.39 
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HR310 
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n 07-5+0-009+0.046 
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U-^0_0.06 -0.11 
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U.iU_o.o6 


p, ,1+0.04 +0.14 
'J-^1-0.04 -0.08 
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<0.15 
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WD 2309+105 
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1 


n 19+0-12 


n n«/1 +0-010+0-052 
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2 
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WD 1210+533 




1 
2 


<0.62 
0.0501°:°^^ 
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p, 07+0.47 +0.57 
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3 


0.071°:?^ 
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2 


o.i4t°:t^ 


p, 10+0.27 +0.31 
U-l°-0.13 -0.14 
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1 


0.28«;?^ 


„ ,.+0.06 +0.11 
U-i4_Q04 -0.06 






(G, Blue) 






2 


0.0l8t°;S°I 


p, ,7+0.01 +0.10 
"--^'-0.01 -0.07 






(Vel) 



^T-wo errors listed: the first are based on the propagation of the column density errors only, while the second include errors in 
the cosmic abundances and the natural range of depletions of carbon and sulfur in the ISM. 

1* Assume LISM tempe rature appropriate for the line of sight based o n multi-ion line w i dths o r the LISM average 
JRedfield fc Linskvll2004bh . For the second component toward 77 UMa, which iRedfield fc Linskvl ||2004b|) estimate a tempera- 
ture of olo'*°° K, we use T = 100 K. 



'^In agreement with projected velocity and spatial distribution IIRedfield &: Linskvll2008h . 
'^Used to calculate the weighted moan value for the LIG, ne(LIG) = 0.12 ± 0.04 cm~^. 
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Fig. la.— Fits of the interstellar C II (1334.5 A) and C IF (1335.6627 A and 1335.7077 A) 
absorption toward 13 nearby stars. The ratio of the column densities are used to estimate 
the electron density for each component. The name of the target star is given above each 
group of plots, and the wavelength (in Angstroms) of each line is provided within each 
plot. Both the 1335.6627 A and 1335.7077 A C 11* lines are shown in the bottom plot. 
Although the 1335.6627 A line is weak, it is evident as the blueward component in some 
spectra (e.g., G191-B2B, 74 Psc A and B, and p Lup). The data are shown in histogram 
form. The thin solid lines are our estimates of the intrinsic stellar flux across the absorption 
feature. The dashed lines are the best-fit individual absorption lines before convolution 
with the instrumental profile. The thick solid line represents the combined absorption fit 
after convolution with the instrumental profile. The spectra are plotted versus heliocentric 
velocity. The parameters for these fits are listed in Table [31 
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Fig. 2a. — Similar to Figure [Ta] but for fits to interstellar S II absorption toward 10 of the 13 
targets which show C 11* absorption. The column density of the optically thin S II profiles 
are used to estimate the column density of C II along the line of sight. All three S II lines 
are fit simultaneously. The parameters for these fits are listed in Table IH 
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Fig. 2b.- 
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Fig. 3. — Similar to Figure [Tal but for fits to interstellar Mg I and Mg II absorption toward 
3 of the 13 targets which show C 11* absorption. The ratio of ionization stages are used 
to estimate the electron density for each absorption component. Both Mg II lines are fit 
simultaneously. The parameters for these fits are listed in Table El 
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Fig. 4. — Estimated C II column density using S II as a proxy [A^(C IIsii)] versus C II reso- 
nance line column density derived from the saturated lines directly [A^(C II)]. The systematic 
errors due to the conversion from A^(S II) (solid red lines) extend beyond the random S II 
fitting errors. Unity (dashed line) bisects the plot window. 
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Fig. 5. — Excited line column density versus electron density (top, black) and versus elec- 
tron density derived from C II column density using S II as a proxy (bottom, blue). The 
systematic errors due to the conversion from A^(S II) (solid red lines) extend beyond the 
random S II fitting errors. Arrows associated with filled symbols indicate errors than extend 
beyond the scope of the plot, whereas arrows associated with open symbols indicate upper 
limits. 
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Fig. 6. — Excited line column density versus resonance line column density (top, black) and 
versus C II column density using S II as a proxy (bottom, blue). The systematic errors 
due to the conversion from A^(S II) (solid red lines) extend beyond the random S II fitting 
errors. Arrows associated with open symbols indicate upper limits. The dotted black lines 
i ndicate lines of equal elec tron density, assuming the LISM average temperature of 6680 K 
(JRedfield fc Linskyll2004bl ). The unweighted mean electron densities are n(,{C U] 



'-'•J^'-'-o.o? 



and ne(C IIsii)= 0.11„o!o5 (thick dashed lines), as calculated from the histograms of He and 
ne(C IIsii) in logarithm (See Figure [7]). 
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Fig. 7. — Histograms of measured electron densities rig (solid, black) and electron densities 
from C II column densities based on S II as a proxy, ne(C IIsn) (dashed, red) in logarithm 
with Gaussian fits. The bin size equals 0.333 dex. The unweighted centroid of the solid 
(black) log-normal distribution is —0.88 log (cm^^) with a dispersion of 0.33 log (cm~^) and 
of the dashed (red) distribution, —0.94 log (cm^^) with a dispersion of 0.26 log (cm~^). The 
shaded (blue) histogram indicates the electron densities of sight lines that are kinematically 
and spatially identified with the LIC All LIC sight lines show a consistent value of rig. 
The top axis gives the estimated pressure P/k = nT, assuming the LISM average value of 
temperature and a simple photoionization relationship between the electron density {rie) and 
the neutral hydrogen density (nni). The axis is not printed for densities in which we expect 
these calculations to fail. For very low densities, the calculate hydrogen density is lower 
than the minimum allowed hydrogen density based on typical observed hydrogen column 
densities and distances to the background star. At very high densities, given a characteristic 
size of a LISM cloud, the resulting high column density (A^hi > 10^^'^), will significantly alter 
the ionization rate due to shielding and therefore contradict the assumption of a constant 
ionization rate. As shown here, the distribution in pressure is simply a function of the 
distribution of the observed electron density. 
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Fig. 8. — Comparison of the temperature dependence in the electron density calculation 
using C ir an d C II (blue) versus Mg I and Mg II (red). This kind of plot was used by 
Gry fc Jenkins! (l200ll ) to put limits on the temperature and electron density of clouds along 
the line of sight toward e CMa. However, an independent measure of the temperature of 
LISM clouds is available from c omparisons of the line widths of ions of different atomic 
mass. iRedfield fc LinskyI (J2004bl ) derive LISM cloud temperatures using this technique and 



the temperature for specific components are shown above by the solid vertical line. The 
gray scale above shows the la, 2o", and 3o" levels of the temperature. By combining the 
information in line widths and ionization abundance ratios, we can more tightly constrain 
both the temperature and electron density of clouds in the LISM. 
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Fig. 9. — left: Measured electron density as a function of the angle from the strongest EUV 
radiation source in the local environment, e CMa. No clear correlation is detected, which 
indicates that the ionization structure of the LISM is not solely determined by this dominant 
source. Instead, it is likely that several of the strong radiation sources in the LISM contribute 
to the distribution of electron density measurements, right: For sight lines in which multiple 
electron densities are measured, we compare the difference in electron density as a function 
of total C II column density. No correlation is apparent, which argues that little shielding 
of ionizing radiation is occuring among the bulk of these clouds. However, there are several 
examples where we can identify dynamically the relative distances of multiple clouds along 
the line of sight, and the more distant cloud (i.e., the cloud closer to the ionizing source) has 
a higher electron density and may shield the nearer cloud, which subsequently has a lower 
measured electron density. 



